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Abstract 

In this paper we set up a rigorous justification for the reinitialization algorithm. Us¬ 
ing the theory of viscosity solutions, we propose a well-posed Hamilton-Jacobi equation 
with a parameter, which is derived from homogenization for a Hamiltonian discontin¬ 
uous in time which appears in the reinitialization. We prove that, as the parameter 
tends to infinity, the solution of the initial value problem converges to a signed distance 
function to the evolving interfaces. A locally uniform convergence is shown when the 
distance function is continuous, whereas a weaker notion of convergence is introduced 
to establish a convergence result to a possibly discontinuous distance function. In terms 
of the geometry of the interfaces, we give a necessary and sufficient condition for the 
continuity of the distance function. We also propose another simpler equation whose 
solution has a gradient bound away from zero. 
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1 Introduction 

Setting of the problem In this paper we establish a rigorous setting for the reinitial¬ 
ization algorithm. In the literature “reinitialization” usually refers to the idea of stopping 
the process of solving an evolution equation regularly in time and changing its solution at 
the stopping time so that we obtain a function which approximates the (signed) distance 
function to the zero level set of the solution. A typical example of such evolution equations 
is 

(1.1) ut = c{x,t)\Vu\, 

where u = u{x,t) is the unknown, Ut = dtU, Vu = ... ,dx„u) and | ■ | stands for the 

standard Euclidean norm in R". The equation fll.ip describes a motion of an interface Pj 
in R” whose normal velocity is equal to c = c{x,t), where at each time the zero level set 
of u{-,t) represents the interface P^. In general, the solution of (ll.l|) does not preserve the 
distance function, and its gradient can get very close to zero. For example, the function 

u{x, t) = 1 — |x| e~^ 

solves the problem 

j Ut = |x| • \ux\ in R X (0,T), 

\ u{x, 0) = 1 — |a:| in R 

in the viscosity sense. For the numerical study of fll.ip several simplifications can be made 
when the solution is or approximates the distance function. One of the reasons is the fact 


2 































that the gradient of the distance function is always 1 and thus bounded away from 0. When 
the gradient degenerates like in the above example, it becomes difficult to compute precisely 
the zero level sets. The reinitialization is used to overcome such an issue. For a more detailed 
discussion on the numerical profits of the reinitialization, see mm- 

Several reinitialization techniques have been introduced in the literature. In this paper 
we focus on the one introduced by Sussman, Smereka and Osher ([Sn])- Their method 
allows to reinitialize fll.ip without explicitly computing the signed distance function with 
the advantage that the level set function of their method approximates the signed distance 
at every time. 

We briefly explain the main idea of the method in [SU]. Consider the corrector equation 

(1.2) (/)* = sign(0)(l - |V0|), 

where sign(-) is the sign function defined as 

{ r 

— if r 7^ 0, 

\r\ 

0 else. 

The solution of this equation asymptotically converges to a steady state | V0| = 1, which is a 
characteristic property of the distance function; see Subsection 13.31 The purpose of the sign 
function in fll.2p is to control the gradient. In the region where (j) is positive, the equation is 
= 1 — |V0|. Thus, the monotonicity of 0 is prescribed by the order of 1 and |V0|. This 
forces I V0| to be close to 1 as time passes. Also, the relation sign(O) = 0 guarantees that the 
initial zero level set is not distorted since 0* = 0 on the zero level. Roughly speaking, the idea 
of [30] is to stop the evolution of (11.11) periodically in time and solve (II.2p till convergence 
to the signed distance function is achieved. This method was first applied in [30| for the 
calculation of the interface of a fluid flow, with the disadvantage that the fluid flow can lose 
mass, because of the accumulation of numerical errors after many periods are completed. 
This problem was later fixed in [25] . 

Up to the authors’ knowledge there is no rigorous setting for the reinitialization process 
described above. In this paper we study an evolution of an interface T^ given as the zero level 
set of the solution u to the initial value problem of the general Hamilton-Jacobi equation 

(1.3) Ut = Hi{x,t,'Vu) 

with a Lipschitz continuous initial datum uq. Here Hi = Hi{x,t,p) is assumed to be contin¬ 
uous, geometric and Lipschitz continuous in x and p. These assumptions are often used in 
the literature to guarantee that fll.31) is well-posed. As a corrector equation we use a slight 
modification of fll.21) . namely 

u 

(1.4) ut = . h(Vu), 

Wo + w 

where £o > 0 is fixed and the function h can be one of the following: 

(1) h{p) = 1 - \p\, 
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(2) The plus part of ([I]), i.e., h{p) = (1 — |p|)+. 


The function /3 (m) = u/\JeQ + is a smoother version of the sign function. Although the 
function h in ([2]) does not preserve the distance function in the sense of [SU] and in a way 
that will be made rigorous later in Theorem 12.21 and Example 14.11 it does however prevent 
the gradient of the solution to approach zero on the zero level set. Moreover, it provides a 
simple monotone scheme for the numerical solution of the problems which we will encounter. 
In fact, our result applies for corrector equations which are more general than fll.4p . but for 
the sake of simplicity we present, in this section, the main idea for this model equation. 

The idea, as in PI. is to solve (11.31) and (ll.dh periodically in time, the first for a period 
of fciAt and the second for k 2 ^t, where fci, /c 2 , At > 0 and one period will be completed at 
a time step of length e = {ki + /c 2 )At. We are thus led to define the following combined 
Hamiltonian 


hfi2(a:,t,r,r,p) := 


if (t - 1) < r < (t - 1) + 


fei At 




■■hiVu) if (i - 1) + ^ < r < i 


for i = 1,..., [j]. Here by [x] we denote the smallest integer which is not smaller than 
X G R. The rescaling of the Hamiltonian Hi in time is required since certain time intervals 
are reserved for the corrector equation. More precisely. Hi is solved in time length kiAt\-^ ~ 


fci 


k\+k2 1 + ■ 


One would expect that solving the two equations infinitely often would force 


the solution of the reinitialization algorithm to converge to the signed distance function to 
T^; we denote it by d. Therefore we are led to study the limit as £ 0 of the solutions of 


( 1 , 6 ) 


ul = Hi2 ( X, t, -, Vn 


H{x, 0 ) = uo{x) 


in X (0,T), 

in R”. 


This is a homogenization problem with the Hamiltonian H 12 being 1-periodic and discontinu¬ 
ous in the fast variable t = t/e. Since the limit above is taken for At —)■ 0 (and consequently 
e —?■ 0), two free parameters still remain, namely ki and k 2 - In fact, we show that the 
solutions of fll.51) converge, as e —)■ 0 and after rescaling, to the solution of 


( 1 . 6 ) 


< = Hi{x, f, Vn^) + ei3{u^)h{\/u^) in R" x (0, T), 
u^{x, 0) = uo{x) in R”. 


Here 6 = k 2 /ki is the ratio of length of the time intervals in which the equations fll.31) and 
(ll.4p are solved. If we solve the corrector equation (11.41) in a larger interval than the one we 
solve the original (11.31) . we can expect the convergence to a steady state. For this reason we 
study the limit as 0 —)■ cx) of the solutions of 01.61) . 

Let us consider the function h in ([T]) as the model case. Roughly speaking, the limit 
9 ^ 00 forces h(Vu^) to be close to 0 except on the zero level of u®, i.e., |Vm®| ~ 1 for large 
value of 6. If we further know that the zero level set of is the same as that of the solution 
of 01.31) and hence is equal to T^ (we call this property a preservation of the zero level set). 
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then we would get a convergence of to the signed distance function d, which is known to 
be a solution of the eikonal equation 

(1.7) |Vd| = 1 

with the homogeneous Dirichlet boundary condition on the zero level. The preservation of 
the zero level set for (11.61) mainly follows from [21]. 

To justify the convergence to d rigorously, the comparison principle for the eikonal equa¬ 
tion fll.7p is used to compare the distance function and a half-relaxed limit of u^, which is 
a weak notion of the limit for a sequence of functions. To do this, we need to know that 
the limit of also preserves the zero level set. This is not clear, despite the fact that 
always preserves the zero level set for every 9 > 0. For the preservation of the zero level 
set by the limit, continuity of the distance function plays an important role. As is known, 
if we £x a time, d{-,t) is a Lipschitz continuous function, but d is not continuous in general 
as a function of {x,t). Indeed, when the interface has an extinction point fPefinition 15.3p . 
the distance function can be discontinuous near this point. For our problem, by construct¬ 
ing suitable barrier functions it turns out that, when d is continuous, the zero level set of 
the half-relaxed limit of is the same as F^. Consequently, we obtain the locally uniform 
convergence of to d; see Theorem 12.11 flilip . 

Concerning the locally uniform convergence, we further consider a condition which guar¬ 
antees the continuity of d. An important property of first order equations is the finite speed 
of propagation (Subsection 15.ip . which allows us to show that the only way the distance 
function can be discontinuous is if points at the zero level extinct instantaneously. More 
precisely, we show that the distance function is continuous at {x, t) if and only if at least one 
of the nearest points of x to Ft is a non extinction point; see Theorem 15.41 ([2]). Therefore, if 
the latter condition is satisfied for every (x, t) G R” X (0,T), then the solutions of fll.61) 
converge locally uniformly to d in x (0,T) fRemark l5.1l) . The converse is also true. 

If the signed distance function d is discontinuous, we cannot expect that the continuous 
solutions of (ll.6p will converge locally uniformly to d. In fact, when d is discontinuous, 
the zero level sets of the half-relaxed limit of are not Ft, and this prevents us to apply 
the comparison principle for fll.7p . We can however show fTheorem l2.1l ([I])) a weaker notion 
of convergence to d; namely a convergence to d from below in time as follows: 

(1.8) lim u^{y,s) = d{x,t) for all {x,t) G R"' x (0,T). 

[y,s,e)^[x,t,oo) 

S<t 

This will be shown by introducing a notion of a half-relaxed limit from below in time and by 
using the fact that d is continuous from below in time. The result fll.Sp also implies a locally 
uniform convergence at any fixed time, that is, converges to d{-,t) locally uniformly 

in R*^ as d —)■ cx); see Theorem 12.11 flill) . 

In a future work we plan to introduce a numerical scheme for (ll.6l) . where no reinitial¬ 
ization will be required. We also plan to study numerically and rigorously a similar method 
for second order equations, including the mean curvature flow, of the form 

Ut = d/i(x, t, Vm, V^m) -|- df3{u)h(yu). 
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Review of the literature In [5], Chopp used a reinitialization algorithm for the mean 
curvature flow. His technique does not utilize a corrector equation, instead, at each stopping 
time he recalculates the signed distance and starts the evolution again with this new initial 
value. In [27] it is mentioned that another way to reinitialize is to compute the signed 
distance at each stopping time using the fast marching method. In [22] the problem of the 
movement of the zero level set which appeared in [SU] is flxed by solving an extra variational 
problem during the iteration of the two equations. Sethian in [27] suggests to use, instead 
of the reinitialization algorithm, the method of extended velocity described in chapter 11. 

In [T2] a new nonlinear equation is introduced for the evolution of open sets with thin 
boundary under a given velocity held. The solution is for every time the signed distance 
function to the boundary of the open set (called an oriented distance function in [12]). Other 
numerical methods for preserving the signed distance are presented in na and ng. 

In [20], the authors use the approximation of the mean curvature flow by the Allen-Cahn 
equation and they prove the convergence of an equation to the signed distance function. See 
also [7] for a related theory developed for anisotropic and crystalline mean curvature flow. 

Summary To sum up, the contributions of this paper are: 

• mathematical justification of the reinitialization procedure, 

• introduction of a new approximate scheme for the distance function of evolving inter¬ 
faces, i.e., solving fll.61) and taking the limit as 0 —)■ cx), 

• formulation of a necessary and sufficient condition for the solution of the scheme fll.bl) 
to converge locally uniformly to the signed distance function, in terms of topological 
changes of interfaces, 

• discovery of a weak notion of a limit which gives the signed distance function even if 
it is discontinuous. 

We also mention that 

• through the rigorous analysis of the reinitialization procedure, we retrieve the correct 
rescaling in time of the equation fll.31) in order to approximate the signed distance 
function, and thus we extend the reinitialization procedure to evolutions with time 
depending velocity fields, 

• lastly, the equation in (11.61) with h satisfying ([2]) or more generally (as we will see later) 
the assumption fl2.1()|) admits a natural numerical scheme with a CFL condition, see 
also [1] or [2B]. We plan to study this last part in a future work. 

Organization of the paper In Subsection 12.11 we state the main results, and in Sub¬ 
section 12.21 we present known results concerning a well-posedness and regularity of viscosity 
solutions. Section [2] consists of main tools which we use in order to prove our main theorems, 
namely the preservation of the zero level set fSubsection 13.ip . construction of barrier func¬ 
tions fSubsection 13.2p and characterization of the distance function via the eikonal equation 
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fSubsection I3.3p . In Section 0] we prove convergence results to the signed distance function 
d. The proof for continuous d and that for discontinuous d will be given separately. Section 
0] is concerned with continuity properties of the distance function, and hnally in Section [5] 
we prove a homogenization result. 

2 Main results 

2.1 Main theorems 

We study the evolution of the zero level set of a function w given by the following problem: 



(2.1a) 

(2,1b) 


Here uq is a possibly unbounded Lipschitz continuous function on R" and its Lipschitz 
constant is denoted by Lq. The function Hi = Hi{x,t,p) : R” x [0,T] x R” —)■ R satishes 

(HI) Hi e ^(R” X [0,T] X R”), 

(H2) Hi{x, t, Xp) = XHi{x, t,p) for all A > 0, x G R”, t G [0, T] and p G R”, 

(H3) There is a positive constant Li such that 

\Hi{x,t,p) - Hi{y,t,p)\ < Li\x - y\ 
for all X, 1 / G R”, t G [0, T] and p G R" with \p\ = 1, 

(H4) There is a positive constant L 2 such that 

\Hi{x,t,p) -Hi{x,t,q)\ < L 2 \p-q\ 
for all X G R”, t G [0, T] and p, g G R"", 

(H5) sup \Hi{x,t,p)\ < 00 . 

(x,t,p)GR’*x[0,T]xR" 

|p|=i 

Remark 2.1. In order to get a more precise estimate for the Lipschitz constant of solutions 
considered in this paper, we will use in Proposition 12.71 instead of flH3j) the following: 

(H3-s) There is a function D G C([0,T]) such that 

\H{x,t,p) -H{y,t,p)\ < D{t)\x-y\ 

for all X, p G R”, t G [0, T] and p G R” with |p| = 1. 

Note that the assumption (lH3-s[) implies the assumption (1H31) . 
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Remark 2.2. In the literature the assumption (IH3j) is usually given as: There are L, L 
positive, such that 

(2.2) \H{x,t,p) - H{y,t,p)\ < L\x - y\\p\ + L\x - y\ 

for all x,y,p G R” and t G [0,T]. However, since the Hamiltonian H is geometric (the 
assumption (lH2D h it turns out that the conditions flH3j) and fl2.2p are equivalent. Indeed, 
it is clear that flH3D implies fl2.2p with L = Li and L = 0. Also, under fl2.2p we can easily 
derive flH3p with Li = L + L. Here let us also show that, in fact, we can take Li = L in flH3D 
when 02.21) holds. Let x^y ^ R"', t G [0,T], p G R"' with IpI = 1 and r > 0. Then we have 

\H{x,t,rp) — iL(|/,t,rp)| < Lr\x — y\ + L\x — y\. 

Dividing both sides by r and using 0H2D . we get 

\H{x,t,p) - H{y,t,p)\ < L\x-y\ + L —— 

r 

If we now take the limit as r —)> +cx), we get 0H3P with Li = L. 

The assumption 0H2p is natural for a geometric evolution problem, while flH3D is used for 
construction of barriers in Subsection 13.2l and for the proof of Lipschitz continuity of solutions 
in Appendix We call the constant L 2 in assumption 0H4p the speed of propagation of the 
zero level set of the solutions. Existence, uniqueness and other properties of the problem 
fl2.ip can be found in Subsection 12.21 Since the zero level set of the solution w of fl2.ll) is the 
main focus of this paper we will use the following notations for t G [0,T): 


(2.3) 

Df := {x E R” 

±w{x, t) > 0}, 

(2.4) 

Tt := {x E R'^ 

w{x, t) = 0} 

and 


D+:= U (D+x{(}), 

D- := U (A'xW) 


te(o,T) 

tG{Q,T) 


In what follows we will always suppose that the evolution associated with w is not empty, 

i.e., 

(2.5) Ti ^ 0 for alH G [0,T). 

Also for D C R" the distance function dist(-, D) ; R" —)■ [0, cxo) is defined as 

dist(a:, D) := inf \x — y\. 
yeo 

Remark 2.3. It is well-known that the distance function dist(-,r2) is Lipschitz continuous. 
Indeed, for x,y E R"" and z G D, we have 

\y — z\ < \x — z\ + \x — y\ and \x — z\ < \y — z\ + \x — y\. 

If we take the infimum for z eVL, the above inequalities become 

|dist(a:, D) — dist(|/, 12)1 < \x — y\. 





















For a function w : R" x [0,T) R, Lip^[tc] stands for the Lipschitz constant of w with 
respect to x, i.e., 


Lip^M 


sup sup 
3;,j/eR" tG[Q,T) 

x^y 


w{x,t) -w{y,t)\ 

\x-y\ 


G [0, cxo]. 


Our hrst result concerns an equation of the form 


(2.6) Vn®)Vn®) inR"x[0,T), 

where 0 > 0 is a parameter, Hi is as in fl2.1ap and 


(2.7) H 2 {r,p) = P{r)h{p). 

The function /3 is assumed to satisfy 

(B) Lip[/3] =: Lp < oo and /? is non-decreasing and bounded in R with /3(0) = 0, /3(r) > 0 
if r > 0, /9(r) < 0 if r < 0, 

where by Lip[/] we denote the Lipschitz constant of a function / ; R” —)■ R. Moreover, 
h : R” —?■ R is such that 

(2.8) there is a modulus Uh such that \h{p) — h{q)\ < Uh{\p — q'|) for all p,q E R"". 

Here a function u : [0, oo) -E [0, oo) is called a modulus if u is non-decreasing and 0 = 
a;(0) = hma;(r). We will also use one of the following assumptions for the function h: 


(2.9a) 

r h{p) > 0 if IpI < 1, 

(2.9b) 

\h(p) < 0 if p > 1 

or 


(2.10a) 

ih{p) > 0 if IpI < 1, 

(2.10b) 

\h(p) = 0 if p >1. 

Examples of these functions 

are 

(2.11) 

Hi{x,t,p) = c(x,f)|p|. 



(2.12) 

H2{%p) - ^ - -h(jp) 


^el + v? 

for 


(2.13) 

Kp) = 1 - bl 
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or 


(2.14) h{p) = (1 - |p|)+, 

where £o > 0, c is Lipschitz continuous with respect to a; G R"' uniformly in time, and for 
a G R we denote by 

a± = max{±a, 0} 

the positive and negative part of a. We see that the function h dehned in 02.131) satishes 
02.9p while 02.14p satishes 02.1UI) . 

For a function w{x, t) dehned in R" x [0, T), we dehne the signed distance function d{x, t), 
from the zero level set of w, as follows: 


(2.15) 


d{x, t) 


dist(a;, Ft) if a; G U F^, 
—dist(a:, F^) if a; G ■ 


Here Df and F^ are dehned in 02.31) and 02.4p . 

For later use we collect our main assumptions in the following list: 


(2.16) 


Mo is Lipschitz continuous in R”, Hi satishes (H1)-(H5), 
H 2 is of the form (2.7), (3 satishes (B), h satishes (2.8). 


For the solution of 02.61) and 02.1bp we have the following main theorem. 

Theorem 2.1 (Convergence of to the signed distance function). Assume 02.161) and 02.91) . 
Let be the solution of 02.61) and 02.1bD . Let d be the signed distance function as in 02.151) . 
Then 


(i) 


lim 

{y,s,e)^{x,t,oo) 

S<t 


u%y,s) = d{x,t) 


for all {x,t) G R” X (0,T), 


(a) u^{-,t) converges to d{-,t) locally uniformly in R"' as 9 ^ +cx) for all t G (0,T), 
(Hi) if in addition d{x,t) is continuous in R"' x (0,T), then 

converges to d locally uniformly in R” x (0,T) as 9 ^ +cx). 


In general, if the signed distance function d is discontinuous, we cannot expect that the 
continuous functions will converge to d locally uniformly. The following example shows 
that the signed distance function can be discontinuous when points of the zero level set 
disappear instantaneously. We will denote by Br{x) the open ball of radius r > 0 centered 
at X. Its closure is Br{x). Also, (•, •) stands for the standard Euclidean inner product. 

Example 2.1 (A single discontinuity). We study fl2.ip with 
(2.17) Hi{x,p) = c{x)\p\, 
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where c G Lip(R"') is bounded and non-negative. Since Hi is written as Hi{x,p) = 
max {c(x)a,p), the viscosity solution w of fl2.ip has a representation formula as a value 

ae_Bi(0) 

function of the associated optimal control problem ([151 Section 10]), which is of the form 

(2.18) w{x,t) = supMo(^“(^))- 

aeA 

Here A := {a : [0,T) —i?i(0), measurable} and : [0, T) —)■ R"' is the solution of the 
state equation 

(X")'(s) = c(X“(s))a(s) in (0,T), X“(0) = x. 

Each element a G ^ is called a control. 

We now consider the case where c{x) = 1. This describes a phenomenon where the 
interface expands at a uniform speed 1. In this case the optimal control forces the state 
to move towards the maximum point of mq in Bt{x), and hence 

(2.19) w{x,t) = niax Uo{y). 

\x y\^t 

Take the initial datum as Uo(x) = max{(l — \x — 2|)_|_, (1 — |x -|- 2|)_|_}. The formula 02.191) 
now implies 

w{x, t) = min{max{(f 1 — |x — 2|)+, (f -|- 1 — |x -|- 2|)+}, 1}. 

To see this we notice that 

max {l-\y± 2|)+ = min{(t + 1 - |x ± 2|)+, 1}, 

\x-y\<t 

then using the formula 02.191) and after changing the order of the maxima, we calculate 

w{x, t) = max{min{(t 1 — |x — 2|)_|_, 1}, min{(t -|- 1 — |x -|- 2|)+, 1} 

= min{max{(f -f 1 — |x — 2|)+, (f -|- 1 — |x -|- 2|)+}, 1}. 

Here we have used the relation max{min{ai, h} min{a 2 , b}} = min{max{ai, 02 }, b} for 
Oi, 02 , 6 G R. We therefore have 

^ Qj ^ I{1^1 - ^ + 3} U {|x| < 1 - f} if f < 1, 

|{|x| > t -|- 3} if f > 1 

and 

. . _ fmax{(f1 - |x - 2|)+, (f-h 1 - |x2|)+} if f < 1, 

’ I (f -|- 3 — |x|)+ if f > 1. 

See Figure m Thus d is discontinuous on £ := {(x, 1) | — 2 < |x| < 2}; more precisely, d is 
not upper semicontinuous but lower semicontinuous on i. 

If h satishes O2.10p we can still estimate from one side the limit with the distance function. 
More precisely we have the following theorem. 
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Figure 1: The graph of d. 


Theorem 2.2. Assume fl2.16p and 02.101) . Then 


d{x,t) < supm{x,t) < +00 

e>o 

dixA) > iniu^ixA) > —oo 
^ ~ e>Q ^ ' 


for all X G Df, 
for all X G Df. 


For the next result we dehne 


( 2 . 20 ) 


Hi 2 {x,t,T,r,p) := 


Hiix,-^,p) if (i - 1) < r < (z - 1) + ^ 


1+^ 


kiAt 


Hiir.p) 


if (i - 1) + 


kiAt 


< T < i 


for ki, k 2 > 0, At > 0, e = (/ci + k 2 )At and z = 1,..., f-]. By dehnition H 12 is 1-periodic in 
r, and in general it is discontinuous in r. In summary, we are led to the following equation: 


(2.21) zzf = ihi 2 t,zz^, Vzz^^ inR"x(0,T). 

Remark 2.4. A solution of the problem 02.2ip . 02.1bP can be constructed by solving 02.211) 
in the intervals [^(z — 1), £(z — 1) + kiAt), [e{i — 1) + fciAt, ez), z = 1,..., [j], iteratively, using 
as initial condition at each interval, the hnal value of the solution dehned in the previous 
interval. We call this solution an iterative solution. 


Let 9 = k 2 /ki. We dehne 

ff{x,t,r,p) = +^^ 2 (ap)) 

and consider the equation 


(2.22) = H{x,t,u\\/u^) inR"x(0,T). 

Theorem 2.3 (Homogenization). Assume 02.16p . Let andu^ be, respectively, the solution 
of 02.22|) . 02.1b[) and the iterative solution of 02.211) . 02.1bp . Then converges to locally 
uniformly in R” X [0,T). 
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Remark 2.5. If we set u^{x, t) = u^{x, (1 + 0)t) in Theorem 12.31 then solves the equation 
fl2.6p and satishes the initial data fl2.1bD . 

Remark 2.6. All of our main theorems have the same assumptions on uq and Hi. For this 
reason, we will assume 02.161) in the rest of the paper except Subsection 13.11 where mq will 
be generalized. For the function h in 02.7p we will differentiate the assumptions 02.9p and 
02.1Up . Finally we will state clearly whether or not the distance function d is continuous. 


2.2 Theorems from the literature 

In this subsection we will present a comparison principle for general equations of the form 


(2.23) Ut = F{x,t,u,'Vu) inR"'x(0,T). 

Let us introduce a notion of viscosity solutions. For this purpose, we first define semicon- 
tinuous envelopes of functions. Let K C R". For a function / : JF —)■ R we denote the upper 
and lower semi continuous envelopes by f* and /* : iF —)■ R U {±C)o} respectively, which are 
as follows: 

f*{z) := hmsup/(|/) = hmsup{/(|/) | y G Bs{z)nK}, 

y^z o—)-0 

f^z) := liminf f{y) = \immi{f{y) \ y G Bs{z) n K}. 

s^o 

Definition 2.4 (Viscosity Solution). We say that u : R"" x [0,T) —)■ R zs a viscosity subso¬ 
lution (resp. a supersolutionj of fl2.23p if u* < +oo(resp. u* > —oo) and if 

<ft < F*{xo, to, (j), V0) {resp. (ft > F^{xo, to, cf, V0)) at Pq = {xo, to) 


whenever 

(u*<(p on BroiPo) ( fu^>(p onBroiPo) 

1 * X . P ] X . P 

[u = (p at Fo V (M* = 0 at Fo 

for (j) G C^(R" X (0,T)), Pq G R" x (0,T) and r > 0 such that Br{Fo) C R"' x (0,T). 

Since we already use the notation D^, we are going to use the symbol for the 
subdifferential respectively for the superdifferential of a function. More precisely for a function 
u : R" X (0, T) —)• R we define a superdifferential J~^u{z, s) of u at {z, s) G R'^ x (0, T) by 


(2.25) 


J^u{z,s) 


|(p,r) G R” X R 


30 G C^(R"' X (0,T)) such that 
(p, r) = {V(j),dt(j)){z,s) and 
max (u — (b) = {u — <i>){z, s) 

R"X(0,T) ^ -ryy , y 


A subdifferential J~u{z,s) is dehned by replacing “max” by “min” in 02.251) . Equivalently, 
we say that a function u : R"" x (0,T) —)■ R is a viscosity subsolution (resp. viscosity 
supersolution) of 02.23P if 


T < F*{z,s, u* {z, s),p) (resp. r > F^{z,s, m * {z, s),p) ) 

for all ( 2 :, s) G R" x (0,T) and (p, r) G J'~^u*{z,s) (resp. (p, r) G J'~u^:{z, s)). 

In order to guarantee the well-posedness of the problem 02.23P and 02.1bD . the following 
assumptions are usually imposed on the function F : R" x 10 , T] X R X R" ^ R. 
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(FI) F e C(R'^ X [0,T] X R X R”), 

(F2) There is an Oq G R such that r i-G- F{x,t,r,p) — a^r is non-increasing on R for all 
{x,t,p) e R” X [0,T] X R”, 

(F3) For R> 0 there is a modulus uji such that 

\F{x,t,r,p) - F{x,t,r,q)\ < ujR{\p-q\) 

for all {x,t,r,p,q) G R” x [0,T] x R x R" x R”, with \p\, \q\ < R, 

(F4) There is a modulus u such that 

\F{x,t,r,p) - F{y,t,r,p)\ < a;(|a: - 2 /|(l + |p|)) 

for all {x,ii,t,r,p) G R” x R"' x [0,T] x R x R"'. 

For the convenience of the reader we will state the comparison principle and sketch its 
proof for the problem (I2.23p and fl2.1bl) . For a detailed proof, see [T^ Theorem 4.1], 

Proposition 2.5 (Comparison pinciple). Assume that F satisfies flFip - flF4H . Let u,v be a 
viscosity subsolution and supersolution respectively of fl2.23l) and assume that they satisfy 

(Al) u*{x, 0) < vfix, 0) for all x G R”, 

(A2) there is a constant K > 0 such that we have on R” x(0,T) 

u{x,t) < K{1 + |a;|), v{x,t) > —K{1 + |a;|), 

(A3) there is a constant K > 0 such that for x,y ^ R" we have 

u*{x, 0) - vfiy, 0) < k\x - y\. 


Then 

u* < u* in R"' X [0,T). 

Proof. Since the comparison principal is more or less classical we will only give a sketch of 
the proof. 

1. We may suppose without loss of generality that u,v are upper, respectively, lower 
semicontinuous. As usual we set 

u(.^t) = e““*M(-,t), v{-,t) = e““*u(-,t), 

where a > Oq and Oq is given by flF2p . Using the notation u, v instead of u, v, we have that 
M, V are sub- and super solutions of the equation 

Mi -t- (a — ao)u = F{x, t, u, Vw), 

where F{x,t,r,p) is non-increasing in r and satishes the conditions flFlD - flF4p . As before we 
denote by F the new F. 
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2. Suppose that 


M= sup {u{x,t) — v{x,t)) > 0. 

(x,pGR"x[0,T) 

We now make the usual doubling of variables trick and define for e,ri, a > 0 
t, y) = u{x, t) - v{y, t) - - a{\x\‘^ + \y\^) 


and 

= sup ^^{x,y,t), 

x,ye'R^,te[ 0 ,T) 

where a = (e, rj, a). As usual we have 0 < < +cx). In order to proceed we need a priori 

bounds on the maximum and to do that we need to be able to control the difference 
u{x,t) — v{y,t) by the modulus \x — y\. One can show (using a doubling of variables trick, 
see for example m Proposition 2.3']) that there is a constant Ct > 0 such that 

u{x, t) - v{y, t) < Ct{ 1 + \x-y\). 

Using this estimate we can show that there is C* > 0 such that 

(y\x\, a\y\ < \JaC. 

The above estimate together with flA3p enables us to find x„,y„ G R” and G (0,T) such 
that ^a{x„,ya,t„) = M„. 

For the term \xa — |/o-p/e, we will need a more refined estimate than the classical one, 
namely we need 

(2.26) lim I limsup = 0. 

£->0 Y ri,a^0 s J 

A proof of a similar estimate can be found in [TOl Proposition 4.4]. 

3. Doubling the variables again in time or using a similar argument as in HD Lemma 2], 
we have for = 2{xa — ya)/^ 

Tf 

— + (a - aQ){u{x^,t„)-v{y^,t^)) 

< F{x^,t^,u{x^,t^),p^ + 2ax^) - F{y^,F,v{y^,F),p^ - 2ay^). 

As in [m Proposition 2.4] there is a 5 > 0 independent of a, such that u{xa, ta)—v{ya, to) > S. 
Using properties of F, flF3ll and flF4F one gets 

Tj 

— + (a - ao)5 < id{\x„ - y„\{l + \pe + 2ax^|)) + UB{2,a{\x„\ + Ip^D), 

where R = o{l jy/e) + o{^/a) as e, a 0. Using 02.261) . we can take the limit first as a, ?7 —)■ 0 
and then as e —)■ 0 to get a contradiction. □ 

Combining Proposition 12.51 with Perron’s method we get the following theorem. 
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Theorem 2.6 (Existence/Uniqueness). Assume 02.161) . Then for all 6 > 0, there exists a 
unique solution u = E C(R"' x [0,T)) of the problem 02.61) and fl2.1bD with 

mR^x[0,T), 

where u'^^{x,t) = Uq{x) + Kt, t) = Uq{x) — Kt, for some K > 0, are a viscosity 

supersolution and subsolution respectively of the same problem. 

Proof. We will only show that UQ±Kt are a subsolution and a supersolution of 02.61) for some 
> 0, large enough depending on since for the rest of the proof we can use a Perron’s 
argument, see for example [23]. We hrst suppose that Uq is smooth. Then by the Lipschitz 
continuity of uq we have |Vmo| < Lq. By assumptions 0H2D and 0H5D . there is a constant 
C > 0 such that |Ri(a;, t, Vmq) | < CLq. Also, since h is continuous and /3 is bounded by 
(!B|) . there is a constant M > 0 such that 

\H 2 {uo,'Vuo)\ = 9\(3{uo)h{'Vuo)\ < 9 max \h{p)\M. 

\p\<Lo 

Finally, if we choose K > 0 such that K > CLq + 9 max \h{p)\M, we get the desired result. 

\p\<Lo 

For the case where uq is not smooth, we use the same argument for elements of the super- 
and sub differential of uq. □ 

The following proposition is proved in Appendix Dehne 

L{t) :=max{Lo,l}e-^“''^'^'''- 

Proposition 2.7 (Lipschitz continuity of solutions). Under the assumptions of Proposition 
\2.A. with Hi satisfying assumption (lH3-s|) instead of (1H31) . the solution u of the problem 
(12.6p and (I2.1bl) satisfies 

\u{x, t) — u{y, t)\ < L{t)\x — y\ for all x,y E R” and t E [0, T). 

Remark 2.7. The Lipschitz continuity of the solution of fl2.ll) will be used in the next 
section to show that the solution of 02.61) gives the same zero level set as 02. ip and that 
there exist barrier functions of independent of 9. There, the Lipschitz constant is allowed 
to depend on the terminal time T. It is well-known that, if Hi is coercive, i.e., 

Hi{x,t,p) -E oo as IpI —)■ cx) uniformly in {x,t), 

then the solution is Lipschitz continuous and its Lipschitz constant does not depend on T. 
See, e.g., [3]. Since such independence of T is not needed for our study, we do not require 
Hi to be coercive in this paper. 

One important property of geometric equations fl2.1al) is the invariance under the change 
of dependent variables. This invariance property as well as the comparison principle play a 
crucial role for the proof of uniqueness of evolutions. 

Theorem 2.8 (Invariance). Let 9 : H ^ H be a nondecreasing and upper semicontinuous 
(resp. lower semicontinuous) function. If w is a viscosity subsolution (resp. supersolution) 
of fl2.1ap . then 9 o w is a viscosity subsolution (resp. supersolution) of fl2.1ap . 
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See [T71 Theorem 4.2.1] for the proof. 


Remark 2.8. As a simple consequence of the invariance property, we see that, when w is a, 
solution of fl2.1ap . the characteristic function on D^{see 02.31) i defined as 



1 if x e Dt, 
0 if x ^ 


is a supersolution of 02.1al) since it is written as x^+(x) = X(o,co)(x) ow. Similarly, X_D+urt(^) 
is a subsolution of 02 .lap . 


It is known that the evolution of the interface {rt}(g(o,T) associated with 02.1al) is inde¬ 
pendent of a choice of the initial data uq. In other words, if the zero levels of initial data are 
the same, then those of the solutions are also the same. See [171 Section 4.2.3 and 4.2.4] for 
the detailed statement and its proof. 


3 Main tools 


3.1 Preservation of the zero level set 


We believe that the preservation of the zero level set is by itself a useful result. For this 
reason we present it in a more general framework than the one we are going to apply it for 
the proof of our main results. 

We study a general equation of the form 


(3.1) 


Ut = t, Vu) + /3(u)G(x, Vu) in R"" x (0, T). 


The function Hi satisfies (1Hll) - flH5l) . For the function /? we assume that ([B]) is true, and for 
G : R” X R" — )■ R we assume 

(G) G satisfies flF3|) . flFip and is bounded from above in R” x R'^. 

Under these assumptions the comparison principle holds for solutions of fl3.ip . Indeed, 
the continuity assumptions on IIi, G and (3 imply that the function F(x, r, p) := Hi{x, t,p) — 
f3{r)G{x,p) satisfies flFlD . flF3D . flF4D while flF2D is fulfilled with 7 = L/^{sup^ny^-^n G). 

To guarantee that solutions of fl3.ll) preserve the original zero level set, two kinds of 
sufficient conditions on G are made in our theorem. One is boundedness of G from below, 
which, unfortunately, excludes the typical case G{x,p) = 1 — |p|. The other condition needs 
only local boundedness of G from below near p = 0 but requires solutions of fl 2 .ip to be 
Lipschitz continuous, which is not true in general if the initial data Uq is just uniformly 
continuous. 

In the first author’s dissertation [7T], Theorem 13.II fill) is established, but we give its proof 
here not only for the reader’s convenience but also in order to show connection with the 
proof of dUp. 
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Theorem 3.1 (Preservation of the zero level set). Let w and u he, respectively, the viscosity 
solution of fl2.ll) and fl3.ip . fl2.1bD with a uniformly continuous Uq. Assume either (p| or (lill) 
below: 

(i) G is bounded from below in R” x R"'. 

(a) G is bounded from below in R” x Bp{0) for some p G (0,1], and Lip^,)^] < oo. 

Then we have P^ = {u{-,t) = 0} and Df = {±u{-,t) > 0} for all t G (0,T), where and 

Pt are defined in fl2.3p and fl2.4p respectively. 

Remark 3.1. Assume that G is independent of x and continuous. Then ([ip is true if G 
satishes fl2.10p . while (jlip is true if G satishes fl2.9l) . 

Proof. Assume that (IIP is true. 

1. Set G* = max{supj^nxRn G, 0} and G* := max{— infRnxRn G, 0}. We dehne 

if w{x, t) > 0, 
if w{x, t) < 0 

if uj(x, t) > 0, 
if w{x, t) < 0 

for {x,t) G R" X [0,T), where Ly is the Lipschitz constant of (3 appearing in (Ep. We claim 
that V* and r* are, respectively, a viscosity supersolution and subsolution of (13.ip . 

2. We shall show that v* is a supersolution. If w is smooth and w{x, t) > 0, we compute 

V* - Hi{x, t, Vv*) = LpG*v* + - Hi{x, t, e^P^^^Vw) 

= LpG*v* + e^P^**{wt - H^{x, t, Vw)} 

> LyG*v* + 0 

> /3(W)G(x, Vo*), 

which implies that v* is a supersolution of (13.11) . In the general case where w is not nec¬ 
essarily smooth, taking an element of the subdifferential of w, we see that v* is a viscosity 
supersolution of fl3.ip . Similar arguments apply to the case when w{x, t) < 0, so that v* is a 
supersolution in {tr > 0} U {tr < 0}. It remains to prove that v* is a supersolution of fl3.ll) 
on {w = 0}. 

Let ( 2 :, s) G R"" X (0,T) be a point such that w{z,s) = 0, and take (p, r) G J~v*{z,s). 
Our goal is to derive 

T > Hi{z, s,p) 

since (3{v*{z,s)) = 0. To do this, we consider a characteristic function g{x,t) = XOt^x)- 
We have v*{z,s) = g{z,s) = 0 and v* < g near (z, s), and thus (p, r) G J~g{z,s). Since 
p is a supersolution of fl2.1ap by Remark 12.81 we have r > Hi{z,s,p), which is the desired 
inequality. Summarizing the above arguments, we conclude that v* is a supersolution of 
fl3.ip . In the same manner we are able to prove that u* is a subsolution of fl3.1l) . 


v*{x,t) := 


e^P^*^w{x, t) 


'w{x, t) 


and 


t>*(a;,f) : = 


e ^P^*^w{x,t) 
e^P^*^w{x, t) 
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3. Since 0) = n*(x, 0) = ^ 0 ( 2 :) for all x G R"", the comparison principle (Proposition 
12.Sp yields 

Vi,{x,t) < u{x,t) < v*{x,t) for all {x,t) G R" x (0,T). 

In particular, we have 

K(-,t) > 0} c > 0} c K(-,t) > 0}. 

Since > 0} = > 0} = by the dehnition of and v*, we conclude that 

= {u{-,t) > 0}. Similarly, we obtain D'j' = {u{-,t) < 0}, and hence Pi = {u{-,t) = 0}. 
Assume that (juj) is true. 

1. Set := max{—infj^„^-g^G, 0} and m := max{Lip,j,[tc], 1}. Instead of u* and v* 

if w{x, t) > 0, 
if w{x, t) <0 

if w{x, t) > 0, 
if w{x, t) < 0. 

Then v* and h* are a viscosity supersolution and subsolution of fl3.ll) respectively. 

2. We shall prove that h* is a subsolution in {tc > 0}. If tc is smooth, we have 

|Vh*| = Vtc| < p, 

which implies that G(a:, Vh*) > —Gi,p. Similarly to Step 2 of (Jlj), we observe 

(h*)i - Ri(x,f, Vh*) = -LpG^pV^ + - Hi{x,t, Vw)} 

< —Lj3Gi,pVi, + 0 

< /3(h*)G(a;, Vh*), 

i.e., h* is a subsolution. The rest of the proof runs as before. □ 

As an immediate consequence of Theorem 13.11 it follows that the evolution which is given 
as the zero level set of the solution of the non-geometric equation fl3.ll) does not depend on 
the choice of its initial data. 

3.2 Barrier functions 

Throughout this subsection we will assume fl2.16p . Thanks to Theorem 13.11 for a general h 
satisfying fl2.8p and either one of the assumptions (12.Qh or (12.101) . the solution of (12.61) and 
fl2.1bp gives the same zero level set as w, i.e., we have P^ = {u^ = 0} and > 0} 

for all t G (0, T). In order to study the behaviour of as 6* —)■ cx) and a relation between P^ 
and the zero level set of the limit of m®, we will construct barrier functions independent of 
6. More precisely we construct an upper barrier /* and a lower barrier /* such that 

/*<«'< r, R = {/* = 0} = {/* = 0}, 

= {/* > 0 } = {/* > 0 }, 

D- = {/* < 0 } = {/* < 0 }. 


dehned in Step 1 of (|I]), we consider the functions 


and 


v*{x, t) : = 


v^{x,t) : = 


{p/m)e 

{p/m)e~^f^^*f‘^w{x, t) 
e^P^*^w{x, t) 
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In this subsection we often use the fact that, if u is a supersolution (resp. subsolution) 
of fl 2 . 6 p in and if F = {u = 0 } and = {±m > 0 }, then is a supersolution 
(resp. subsolution) of (l2.b|) in R"' x (0,T). This follows from Remark 12.81 Indeed, if 
(p, r) G J~u+{z, s) and u+{z, s) = 0, then we have (p, r) G J~Xd+{.z, s) and this yields the 
desired viscosity inequality since the characteristic function is a supersolution of fl 2 . 1 al) by 
Remark 12.81 The proof for a subsolution is similar. 

We first show that the solutions are monotone with respect to 9 when h is nonnegative. 
This gives a lower barrier in and an upper barrier in D~ in the case of (| 2 . 1 (J|) . 

Proposition 3.2 (Monotonicity). Assume that h > 0. Let 0 < 9i < $2 and and be, 
respectively, the viscosity solution of fl2.6p with 9 = 9i and 6 * 2 . Then 

(3.2) u^^{x,t) < u^^{x,t) for all X E Df, 

(3.3) u^^{x,t) > u^^{x,t) for all X E Df. 

Proof. In D'*' we observe 


= Ri(x, Vm^i) +0i/3(M®i)h(VM^i) 

< Hi{x, Vm^i) + 92(5{u^^)h{yu^^) 

since h is nonnegative. This implies that is a subsolution of fl2.6p with 9 = 92 in D~^. 
Applying the comparison principle to a subsolution (m®^)+ and a supersolution (m®^)+ of 02.61) 
with 9 = 92, we conclude in D^. By the same argument we see that in 

D-. □ 

We show that solutions of 02 . lap with small Lipschitz constants give rise to lower barrier 
functions in and upper barrier functions in D~. 

Proposition 3.3. Assume that h{p) > 0 if \p\ < 1. Let w be the solution of 02.11) . Then 
the viscosity solution of 02. 6 p and 02.1bD satisfies 

(3.4) u^{x,t) > sw{x,t) for all X E Df , 

(3.5) u^{x,t) < ew{x,t) for all x E Df, 

where e := min{l/Lip,j,[t(;], 1}. 

Proof. Set w := ew. Since |Vti)| = £|Vtc| < 1, by the assumption of h we observe 

Wt = Hi{x, t, Vw) < Hi{x, t, Vw) + l3{w)h{'\/w) 

if tc > 0. In other words, tc is a subsolution of 02.61) in {w > 0}. Applying the comparison 
principle to a subsolution (tc)+ and a supersolution (m®)+ of 02.61) . we obtain 03.4p . The 
estimate 03.5p is shown in a similar way. □ 

It remains to construct an upper barrier in and a lower barrier in D~. In both the 
cases 02 .9p and O 2 . 10 p . the solutions are dominated by the signed distance function d with 
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large coefficient. In the proof of Proposition 13.41 below, we use the fact that h is a viscosity 
super solution of 

(3.6) dt = Hi{x — d'Vd,t,'Vd) in {d > 0}. 

This assertion is more or less known (see, e.g., [HI Proof of Theorem 2.2, Step 1-3]), but we 
give its proof in Remark 13.21 for the reader’s convenience. 

Proposition 3.4. Assume that h{p) < 0 if \p\ > 1. Then the viscosity solution of 02.61) 
and 02.1bD satisfies 

(3.7) u^{x,t) < le^^^d{x,t) 

(3.8) u^{x,t) > le^^^d{x,t) 

where I := max{Lo, 1} md Li is the constant in 0H3I) . 

Proof. Dehne d{x,t) := le^^*d{x,t). If d is smooth, then 

dt — Hi{x, t, Vd) = ILie^^^d + le^^^dt — Hi{x, t, le^^^Vd) 

= le^^^{Ltd + dt - Hi{x,t,Vd)}. 

We next apply the fact that d is a supersolution of 03.6p to estimate 

dt — Hi{x, t, Vd) > le^^*{Lid + Hi{x — dVd, f, Vd) — idi(x, f, Vd)} 

> le^^^{Lid - Li\dVd\\Vd\} 

if d > 0. Noting that |Vd| = 1, we have 

dt — Hi{x, t, Vd) > le^^^{Lid — Lid} = 0. 

Since |Vd| = /e^i*|Vd| > 1, we now have h(Vd) < 0 by assumption. This implies that d is a 
supersolution of 02.6p in {d > 0}. Even if d is not smooth, the same arguments above work 
in the viscosity sense. 

Finally, since Uq < ld+{-, 0) in R”, applying the comparison principle to a subsolution 
and a supersolution (d)+ of 02.61) . we conclude 03.7p . The proof of 03.8p is similar. □ 

Remark 3.2. We shall explain why d is a supersolution of 03.61) . We hrst note that d is lower 
semicontinuous in D'^ fTheorem 15.41 OTDi. Let (xo,to) G R" x (0,T) be a point satisfying 
d(xo,to) > 0 and take any (p,r) G J'~d{xo,to)- We choose a smooth function 0 G such 
that (p, r) = (V0, 0i)(a:o,to) and 

min (d — 0) = (d — 0 ) (xq, to) = 0. 

R’*x(o,r) 

Set do := d(xo,to)- Since p G j7“(d|t=tg)(xo), it follows that the closest point of Ti^ to xo 
is unique and that this point is given by po := ^^o — dop G Ftp; for the proof, refer to [21 
Proposition II.2.14] or [HI Corollary 3.4.5 (i), (ii)]. We also remark that |p| = 1. 


for all X G Df, 
for all X G Df, 
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Define '^(x, t) := 0(a: + dop, t) — do. We now assert 


(3.9) 


min 

R’*x(0,T) 


{d+ - V') 


{d+ - 'ip){yo,to). 


Since {d+—'ip){yo, to) = 0 and d+ > 0, we only need to show {-0 > 0} C {d > 0}. Take a point 
{x,t) G R"" X (0,T) such that '^(x, t) > 0. We then have d{x + doP,t) > (j){x + doP,t) > do- 
Using the Lipschitz continuity of d, we compute 


d{x, t) > d{x + doP, t) — do\p\ > do — do = 0. 

Thus fl3.9p is proved. Let g{x,t) = x^+(x). Then the relation fl3.9p implies that (p, r) G 
xJ~ g{yo,to), where we applied (VV', V't)(2/o, ^o) = (V(/), 0t)(xo, to) = {Vi't)- Since the charac¬ 
teristic function p is a supersolution of fl2.1ap (see Remark I2.8p . we have 

T > Hi{yo,to,p) = Hi{xo - dop,to,p), 

which is the inequality we need in order to conclude that d is a supersolution of (l3.6l) . 

Remark 3.3. Another way of proving Proposition 13.41 is using the Lipschitz continuity of 
solutions of (l2.6p and fl2.1bp from Proposition 12.71 Using assumption flH3p instead of flH3-sD 
in Proposition 12.71 the Lipschitz estimate for reads as follows: 


(3.10) \u^{x,t) — u^{y,t)\ < le^^^\x — y\ for all x,p G R", f G [0,T), 

where I = max{Lo, 1} (-^o = Lip[Mo]). If we take the inhmum for all y G Pf in 03.101) we get 
— /e'^^*dist(x, Pt) < u^{x,t) < /e^^*dist(x, Pi) for all {x,t) G R" x [0,T), 
which implies the relations 03.7p and 03.81) . 


3.3 Comparison principle for eikonal equations 

We investigate uniqueness of solutions of the eikonal equation |Vm| = 1 in a possibly 
unbounded set. To establish a convergence to the signed distance function, we show in 
the next section that the limit of the solutions solves the eikonal equation. Since the 
distance function is a solution of the eikonal equation, the uniqueness result presented below 
guarantees that the limit is the distance function. 

We consider the eikonal equation 

(3.11) |Vm| = 1 in D 
with the boundary condition 

(3.12) u = 0 on 512. 

Here 12 C R"^ is a possibly unbounded open set. We denote by d^ the distance function to 
512, i.e., dn{x) := dist(x, 512). It is well known that dn is a viscosity solution of (Id.llh : see. 
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e.g., [21 Corollary 11.2.16] or P, Corollary 3.4.5 (i), (ii) or Remark 5.6.1]. In other words, 
the problem 03.111) with 03.121) admits at least one viscosity solution. Comparison principle 
(and hence uniqueness) of viscosity solutions of 03.11|) and 03.121) is established in [21] when 
is bounded. If is not bounded, the uniqueness of solutions does not hold in general; for 
instance, when = (0, cx)) C R, all of the following functions are solutions: 

dn{x) = X, —dn{x) = —x, Ua{x) = min{x, a — x} (a > 0). 

However, even if H is not bounded, it turns out that nonnegative solutions of 03.111) and 
03.12P are unique and equal to 

Lemma 3.5. Let m H —)■ R. 

(1) If u is a viscosity subsolution of 03.11|) and u* <0 on dfl, then u* < d^ in fl. 

(2) If u is a viscosity supersolution of 03.lip and u > 0 in then < m* in H. 

Proof. (1) It is known that every subsolution of 03.lip is Lipschitz continuous with Lipschitz 
constant less than or equal to one, that is, \u*{x) — u*{y)\ < \x — y\ for all x,y E Q. (For the 
proof see, e.g., [IHl Lemma 5.6] or [251 Proof of Proposition 2.1, Step 1].) This yields the 
inequality u* < d^. 

(2) We consider abounded set H/j := QnB^^O) with R > 0. Dehne d/j(x) := dist(a;, 

We hrst note that > 0 on H since m > 0 in H, and that > 0 = d^j on dQ^. Thus, by 
the comparison principle in bounded sets, we see d/j < in Qr. Finally, sending R —)■ cx), 
we conclude d^ < m* in fl. □ 

4 Convergence results 

Throughout this section we assume (12.161) . We will hrst prove Theorem 12.11 (lilip . it will 
then be easier for the reader to understand the proof of Theorem 12.11 (Ej), (Ell). 

We introduce a notion of the half-relaxed limits ([IDl Section 6]), which are weak limits 
of a sequence of functions and will be used in the proof of the convergence to the distance 
function. We dehne an upper half-relaxed limit u = lim supg_,.^ and a lower half-relaxed 
limit u = lim inf *6i_j.oo as 

u{x,t) := limsup u^{y,s) 

{y,s,0)^{x,t,oo) 

= hmsup{M^(|/, s) I lx — |/| < d, |f — s| < d, d > 1/d}, 

(5—>-0 

u{x,t) := liminf u^{y,s) 

(y,s,d)^(x,t,oo) 

= lim inf{M^(|/, s) | |x — |/| < d, |f — s| < d, d > 1/d}. 

(5—>-0 

Thanks to the existence of barrier functions shown in Section 13.21 we see that, in both the 
cases 02.91) and O2.10p . —oo < u < oo and —oo < u< oo. 

The following proposition is true in the general case where the distance function is not 
necessarily continuous. 
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Proposition 4.1 (The zero level set of the relaxed limits). Assume either 02.91) or 02.1UI) . 
Then 


(4.1) 

{m > 0} = D+, 

{m = 0} C F, 

{m < 0} D F 

and 




(4.2) 

{m > 0} D D^, 

{m = 0} C F, 

{m < 0} = D~. 


Proof. We only show 04.ip since a proof of 04.2p is similar. Let v ;= liminf*e^oo(w®) + . Then 
it is easily seen that v = («)+. From the estimates 0d.4p and 0d.7p of by barrier fnnctions 
we derive 

ew+ < {u^)+ < Ld+ 

for some L > 0. Taking the lower half-relaxed limit, we obtain 


ew^ < V < L((i_|_)* < Ld+. 


Since > 0} = {d+ > 0} = the above ineqnalities imply {n > 0} = D'^, and hence 
{m > 0} = D~^. We similarly have 


—Ld^ < —(m®)_ < —ew-. 

In this case, however, taking the lower half-relaxed limit yields only {m < 0} D {w < 0} = 
D~ becanse —d- is npper semicontinnous. The inclusion {m = 0} C T is now clear. □ 

4.1 Convergence results for continuous distance function 

The following general properties of the relaxed limits will be used to prove the convergence 
of u^-. 


• Assume that each is a subsolution (resp. supersolution) of the equation Fg = 0. If 
Fg converges to some F locally uniformly and u < oo (resp. u > —oo), then u is a 
subsolution (resp. n is a supersolution) of F = 0. 

• \iu = u=\u and —cx) < u < oo, then converges to u locally uniformly as 6* —)■ cx). 

See [TUI Lemma 6.1, Remark 6.4] for the proofs. 

Assume that d is continuous in R” x (0,T), in particular, we can now use the additional 
upper-semicontinuity property of d+ and d-. Proceeding in a similar way as in Proposition 
14.11 we can show 

(4.3) P = {« = 0} = {m = 0}, = {±u > 0} = {±M > 0}. 

Proof of Theorem \2.1\ (Inip . 1. For 0 > 0 we dehne 

Fg{x,t,r,p,T) := ^{r - Hi{x,t,p)} - P{r)h{p). 
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Then is a viscosity solution of the equation Fg{x, t, u, Vm, Ut) = 0 in R” x (0, T). Since Fg 
converges to —(5{r)h{p) locally uniformly as 6^ —)■ oo, it follows that u and u are, respectively, 
a viscosity subsolution and a viscosity supersolution of —(3{u)h{'Vu) = 0 in R" X (0,T). 

Recall that h satisfies fl2.9p . Since f3(u) >0 in D~^ and P(u) < 0 in D~ by fl4.3p . we see 
that M is a subsolution of 

(4.4) \Vu{x,t)\ = l in 
and 

(4.5) —\Vu{x,t)\ = —l in D~ 

as a function of {x,t). (Note that these two equations are different in the viscosity sense.) 
Similarly, m is a supersolution of both (14.4p and (14.51) . Thus, for each fixed to G (0, T), u\t=to 
and u\t=to are, respectively, a subsolution and a supersolution of 


(4.6) |Vm(x)| = 1 inD+ 

as a function of x. (See Remark 14.11 for the details.) By Lemma [3.51 we obtain 


u\t=to < d{-,to) < u\t=to in At 


and hence 


d = u = u in . 


This implies that ^ d locally uniformly in . For D~ we notice that if u{-,to) is a 
subsolution of —|Vm| = —1 then —M(-,to) is a supersolution of |Vm| = 1, hence a comparison 
with —d this time gives the desired result. □ 


Remark 4.1. We claim that, if m = u{x,t) is a subsolution of 04.41) . then u\t=to is a subso¬ 
lution of 04.6p for a fixed to G (0, T). To show this, we take a test function 0 G C^(R"') such 
that maxRn(M|i= 4 g — 0) = u{xo, to) — 0(xo) for xq G We may assume that this is a strict 
maximum. Next define 'ijjM{x,t) := 0(x) -|- M{t — fo)^- We then have 


(lim inf*0M ) (x, t) 
V M—¥oo ' 


0(a:) if t = to, 
oo if t 7 ^ to, 


so that u — (liminf* 0 m) has a strict maximum over R"' x (0,T) at (a:o,to). By [T71 Lemma 
2.2.5] there exist sequences C (0, oo) and {{xn,tn)}^=i C R"" x (0,T) such that 

Mn — )■ oo, {xn, tn) — t (a^o, to) as u — >• OO and u — 'ipM^ has a local maximum at {xn, tn)- Since 
M is a subsolution of 04.41) . we have 


1 > |V0M„(a:^n,i:n)| = |V0(a:„)|. 

Sending n ^ oo implies |V0(xo)| < 1; namely, u\t=to is a subsolution of 04.6p . 
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4.2 Convergence results for general distance functions 

For the case where the distance d is not necessarily continuous we can only compare the 
half-relaxed limits with the distance function in certain domains due to the fact that only 
the inclusions in Proposition 14.11 are true. 

Lemma 4.2 (Comparison with the distance). Assume that either fl2.9p or fl2.10p hold. Then 

(0 

(4.7) 

(4.8) 

(2) For every t G (0, T), 

(4.9) u{-,t)=0 on dDf 

(4.10) M(-,t)=0 on dDT■ 

Proof. We give proofs of 04.71) and 04.9p since 04.81) and 04.101) can be shown in similar ways. 

([T]) In the same manner as in the proof of Theorem 12.11 fllnl) . it follows that u{-,t) is a 
viscosity supersolution of 04.Op in Df. Since u{-,t) > 0 in Df by 04.ip . the comparison 
principle iLemma 13.51 fl^l implies that d{-,t) < lL{-,t) in Df. 

(12]) By 03.71) and 03.8p we have 

(4.11) -Ld_ <u^ < Ld+, 

where L > 0 is a constant. Taking the lower half-relaxed limit at {x,t), we obtain 

(4.12) —L{d-)*{x,t) < u{x,t) < L{d+)^,{x,t) < Ld^{x,t)- 

Let X G dDf. Then the right-hand side of 04.12p is 0 since x G T^. We next study the limit 
of (i_(|/, s) = dist(|/, Df U Ts) on the left-hand side. Since x G dDf C Df C mt{Df U T*), 
it is not an extinction point (Dehnition 15.31) by Proposition 15.51 Therefore Theorem 15.41 (15p 
ensures that d_ is continuous at {x,t)- This implies that the left-hand side of 04.121) is 0, 
and hence the conclusion follows. □ 

Proof of Theorem \2.‘A Since (I2.1()|) holds, the monotonicity of (Proposition 13.21) yields 
the following representations: 

u{x,t) ={x,t) for X G m(x, f) = inf m®(x, f) ioi x & Df. 

e>o 

These relations and Lemma 14.21 (fip conclude the proof. □ 

For the equation with h satisfying 02.101) . Theorem 12.21 guarantees only the one side 
inequality between the supremum of and the signed distance function d. However, as the 
next example shows, the opposite inequality is not true in general even if the initial datum 
is smaller than the distance function. 


d < u in 
u < d in D~. 
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Example 4.1. Let us consider fl2.ip for the Hamiltonian of the form 02.17p with c(x) = 
(1 — |a;|)+ + 1. We take the initial datum uq as Uo(x) = (1 — |x|)+. The unique viscosity 
solution w of this initial value problem is given as the value function 02.18p . In this case the 
optimal control is the one that leads to a straight trajectory with the maximal speed before 
it comes to the origin and stays there after that moment. Thus direct calculations yield the 
following simplified representation of w: 


w(x, t) 


'1 

(2 — \x\)e^ — 1 
gt-hl+i _ \ 

0 

\ 


if |a:| < 2(1 — e *), 
if 2(1 -e-*) < |a:| < 1 , 
if 1 < |x| < t + 1 , 
if t + 1 < |x| 


for t < log 2 , 


and 


fl if |a:| < t + 1 — log 2 , 

w{x^ ^ I — 1 if t + 1 — log 2 < |x| < t + 1 , for t > log 2 . 

[o if t + 1 < |a:| 


See Figure [21 In particular, we have w{x,t) = 1 if |a:| = t + 1 — log 2 > 1. Also, {tc = 0} = 
{|x| > t + 1 } and the signed distance function d to the interface is d{x,t) = (t + 1 — |x|)+. 
We thus have d{x, t) = log 2 if |x| = t + 1 — log 2, and so 


(4.13) d{x,t) = \og2 < 1 = w{x,t) if |a;| = t + 1 — log2 > 1. 


Since the solution w is non-negative, it is a viscosity subsolution of fl 2 . 6 p with h > 0 for 
every 6 > 0. Accordingly, w < hy the comparison principle. From 04.131) it follows that 

d{x, t) = log 2 < 1 < u^{x,t) if |a;| = t -t -1 — log 2 > 1 , 

which implies that the inequality d > supg^o does not hold on the whole space. 

We also remark that, for 7 G (log2,l), the inequality d{x,t) < yw{x,t) holds if |a:| = 
t -|- 1 — log 2 > 1 and that yw is a solution of (12.lap with the initial datum yuQ. From this 
we see that can be greater than d at some point even if we take an initial datum which 
is strictly less than d{x, 0 ) in {d{-, 0 ) > 0 }. 



Figure 2: The graph of w when t > log 2. 
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In the rest of this subsection we will assume that h satishes the assumption fl2.9p . We 
now introduce several notions of half-relaxed limits. Let {x,t) G R" x (0,T). We define an 
upper and a lower half-relaxed limit from below in time by, respectively, 

vf{x,t):= limsup u^{y,s), yf{x,t) := liminf u^{y,s). 

s<t 

An upper and a lower half-relaxed limit at a fixed time are, respectively, given as 
u\t{x) := limsup u^{y,t), u\t{x) := liminf u^{y,t). 

iy,e)^(x,oo) - (y,e)^ix,oo) 


By definitions we have 


(4.14) u{x,t) < u'{x,t) < u\t{x) < u\t{x) < u'{x,t) < u{x,t) 

for all {x,t) E R"" x (0,T). 

The next proposition is a crucial step in proving a convergence to the signed distance 
function in a weak sense. 


Proposition 4.3. Assume either 02.91) or 02.1UI) . Then the functions u'{-,t) andy!{-,t) are, 
respectively, a viscosity subsolution of 04.4p in Df and a viscosity supersolution of 04.51) in 
Df for every t E (0,T). 

Proof. Fix t E (0,T) and let us prove that u'{-,t) is a viscosity subsolution of 04.4p in . 

1. We hrst introduce an upper half-relaxed limit of in R” x (0, f]. For {x, t) E R"" x (0, t\ 
we define 

v{x,t) := limsup u^{y,s), 

s<t 

which is an upper semi continuous function on R” x (0,t]. By definition we have 


u{x, t) a t < t 
u'{x, t) a t = t. 


2. Take z E and E C'^(R"') such that u'{-, t) — ip attains a maximum at z over R"". 
As usual we may assume that this is a strict maximum, and note that, by 04.71) and 04.14p . 


(4.15) 

We now define fP{x,t) := 


0 < d{z, t) < u{z, t) < u'{z,t). 
ip{x) — \/6{t — i) and 

-|-oo ii t <t, 


(p{x,t) := 


ip{x) ii t = i. 


Then v — (p attains its strict maximum at (z, t) over R” X (0, f], and — (p^ ^ v — (p in the 
sense of the upper half-relaxed limit on R" x (0,f]. Thus, by [T71 Lemma 2.2.5] there exist 
















sequences C (0, oo) and C R” x (0,f| such that 6j —?■ cx), {xj,tj) {z,t) 

and tj) {v — 4>){z, t) as j —)■ cx). 

We now claim 

(4.16) u\z,i) = lim u^^(xj,tj). 

J-fOO 

Observe 

u^'{xj,tj) = {{u^^ - 4)^^){xj,tj) - {v -4)){z,i)} + 4)^^{xj,tj) + (y-4)){z,t) 

= {{u^^ - (j)^^){xj,tj) - {v-(p){z,t)} + {'ilj{xj) - 'ijj{z)} + if{z,t) - - t) 

> {{u^^ - (p^^){xj,tj) - {v - (p){z,t)} + {ip{xj) -ip{z)} + u'{z,t}. 

This implies liminf > u'{z,t). The opposite relation limsupu®-’(xj, tj) < vf{z,t) 

j^oo 

follows from the dehnition of u', and therefore (14.161) is proved. 

3. Since is a viscosity solution of (12.61) in R” x (0,t) and since the viscosity property 
is extended up to the terminal time t = t (0 Section 7]), we have 

(th'{xj,tj) < Ri(xj,tj, V/^(xj,tj)) + 6»^(M®^ (xj,tj))h(V/^(xj,tj)). 

By the dehnition of 0®, this is equivalent to 

—\/9< Hi{xj,tj,'V'il^{xj)) + 9f3{u^^{xj,tj))h{'V'ijj{xj)). 

Dividing both the sides by 6 and sending 6 ^ oo, we obtain 

0 < l3{u{z,t))h{V'ijj{z)), 

where we have used fl4.16p . Since (3{u'(z,t)) > 0 by (I4.15p . using the assumption on h, we 
conclude that |V' 0 ( 2 ;)| <1. □ 

As a consequence of Proposition 14.31 we obtain 

Theorem 4.4. Assume fl2.9l) . Then the following hold. 

(1) u{-, t) = d{-, t) on Df and u{-, t) = d{-, t) on Df for every t E (0, T). 

(2) u' = yf = d in R” x (0,T), i.e., 

lim u^{y,s) = d{x,t) for all {x,t) E R” x (0,T). 

{y,s,e)^{x,t,oo) 

S<t 


(3) u\t = u\t = d{-,t) in R" for every t E (0,T), i.e., u^{-,t) converges to d{-,t) locally 
uniformly in R"" for every t E (0,T). 
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Proof. 1. We first note that fl4.1ip yields 


(4.17) u'= v'= 0 on r, 

(4.18) u\t = u\t = 0 on r. 

Indeed, for {x, t) G F, taking the npper and lower half-relaxed limit from below in time in 
we see that u'{x, t) = u'{x, t) = 0 since d is continuous from below in time by Theorem 
15.41 (121). Similarly, 04.181) follows from the continuity of Thus (ED and (ED were proved 

on r. The equalities in ([ID on dDf or dDf are consequences of Lemma [4.21 (Ep . 

2. It remains to prove ([ID-(|2D in Df and Df. Recall that and are, respec¬ 

tively, a viscosity subsolution of 04.4p in Df and a viscosity supersolution of 04.5p in Df by 
Proposition 14.31 Since 04.171) holds, the comparison result (Lemma 13.51 (ITDl implies that 

(4.19) < d{-,t) inDf, 

(4.20) d(-, t) <«'(•, t) in Up. 


Combining 04.71) . 04.14p and 04.191) . we obtain 


0 < d{-, t) = «(•, t) = «'(•, t) = u\t = u\t = u'{- 


in Df. 


In the same manner, we see 




0 > d{-,t) = = u\t = u\t = = «(• 

■,t) 

in Df. 


The two relations above conclude the proof. 



□ 


This concludes the proof of Theorem 12.11 (ED and fllll) . 


5 Continuity of distance functions 

Throughout this section we study only non-negative distance functions. Namely, we 
assume D~ = 0 so that d{x, t) = dist(x, R) > 0 for all (x, t) G R” x (0, T). Also, we simply 
write Dt = Df and D = . In the general case where d can take negative values, we 

decompose d a,s d = d+ — d_ and apply the following results to d^ and d_. 

5.1 Finite Propagation 

In order to study the continuity of distance functions, we hrst prepare a property of hnite 
propagation for the Hamilton-Jacobi equation fl2.ll) . For this property, the assumption flH4p . 
the Lipschitz continuity of Hi in p plays an important role, though we omit the details in 
this paper. 

Let (x, t) G R"' X (0, T) and r > 0. We dehne a cone as 

CUt)--= U i?.-r(x)x|t+f|. 

0<r<r ^2 J 
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Theorem 5.1 (Local comparison principle). Let {x,t) G R" x (0,T), r > 0 and set C : = 
Uu,v G C{C) are, respectively, a viscosity sub- and supersolution of fl2.1al) in C and 
w(-, t) < v{-, t) in Br{x), then u < v in C. 

See [21 Theorem III.3.12, (Exercise 3.5)] or in [H Theorem 5.3] for the proof. As a 
consequence of Theorem 15.11 we obtain 

Proposition 5.2 (Finite propagation). Let {x,t) G R"" x (0, T) and r > 0. 

(1) If Br{x) C Dt, thenC'^^ ^-^ C D. 

(2) If Br{x) C Ft, then C T. 

Proof. Let w be the solution of fl2.ll) . 

o Set a := mi n tn > 0, and dehne u{x,t) := a, which is a constant function 

Br{x)x{t} 

satisfying u{-,t) < w{-,t) in Br{x). Moreover, u is a solution of (12.lap by the geometricity 
of Hi. Therefore Theorem 15.11 implies that u < w in The positivity of u implies the 

conclusion. 

([2]) The proof is similar to ([T]). We compare w with u{x,t) := 0 both from above and 
from below to conclude that 0 = M<tn<M = 0in ^. □ 

5.2 Continuity properties 

We hrst introduce a notion of extinction points. 

Definition 5.3 (Extinction point). Let x G T^. We say that x is an extinction point if there 
exist 8,5 > 0 such that B^{x) x {t,t + 5] C D. 

For example the point 0 G Ti in Example 12.11 is an extinction point. We remark that 
X G Ti is non-extinction point if and only if there exists a sequence {{xj,tj)}°Li such that 
{xj,tj) — )■ (x,t) as i oo, Xj G Tt^. and tj > t for all j. Dehne Et C R” as the set of 
all extinction points at time t G (0,T) and Nt{x) as the set of all the nearest points from 
X G R” to Tt, i.e., 

Nt{x) := {z G Tt I d{x,t) = \x — z\}. 

Note that we always have Nt{x) 7 ^ 0 by fl2.5p . 

Theorem 5.4 (Continuity properties of the distance function). (1) d is lower semicontin- 
uous in R"' X (0, T). 

(2) d is continuous from below in time, i.e., 

d{x,t) = lim d{y,s) /or aH (x, t) G R" x (0, T). 

(y,s)^{x,t) 

s<t 

(3) Let (x, t) G R” X (0,T). Then d is continuous at (x, t) if and only if Nt{x) \Et (/). 
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Proof. ([T]) The proof can be found in [TTl Proposition 2.1], 

([2]) 1. Suppose by contradiction that d is not continuous at {x,t) from below in time. 
Since d is lower semicontinuous by ([I]), we would have a sequence {{xj,tj)}fLi such that 
{x, t) as j —)■ oo, tj < t and 

lim d{xj,tj) > d(x,t). 

j^ao 

Set a := d{xj,tj) — d{x,t)}/4: > 0. Without loss of generality we may assume that 

d{xj, tj) — d{x, t) > 3a and \xj — x\ < a for all j > 1. 

2. Take any z e Nt{x). We claim that 


(5.1) d{y,tj) > a for all y G Ba{z) and j > 1. 

Since d{-,tj) is a Lipschitz continuous function with the Lipschitz constant 1, we calculate 


d{y,tj) > d{xj,tj) - \xj - y\ 

> {d{x, t) + Sa} — {\xj — x\ + \x — z\ + \z — y\) 

> {d{x, t) + 3a} — (a + |x — .^1 + a) 

> a, 

which yields fl5.ll) . By fl5.ip we have Ba{z) x {tj} C D. Thus Proposition 15.21 (fT|) implies 
that 


(5.2) C D. 

Since tj f t a,s j ^ oo, we have {z,t) G C^zt ) j large, and therefore z E Dt hj 05.2|) . 
However, this contradicts the fact that z EVt- 

([3]) 1. We hrst assume that d is continuous at {x,t). Take any sequence {{xj,tj)}jPi 
such that {xj,tj) —)■ {x,t) as j —)■ cx) and tj > t. By continuity we have d{xj,tj) —)■ d{x,t) 
as j —)■ cx). We now take Zj G Nt.[xj) for each j. Then {zj} is bounded. Indeed, since 
\zj\ < |x| + \x — Xj\ + \xj — Zj\ and \x — Xj\ —)■ 0, \xj — Zj\ = d{xj, tj) —)• d{x, t) as j —)■ oo, we 
see that {zj} is bounded. From this Zj subsequently converges to some z as j —)■ oo, where 
we use again the index j. It is easy to see that z ETf. 

Let us show z G Nt{x) \ Et. Taking the limit in d{xj,tj) = \xj — Zj\, we obtain d{x,t) = 
\x — z\, which implies that z G Nt{x). Also, since Zj E Pt- and tj f t as j ^ oo, it follows 
that z is not an extinction point, and hence we conclude that Nt{x) \ Et 

2. We next assume that d is not continuous at {x,t). By ([T]) and ([2]) we have some 
sequence {{xj,tj)}'jLi such that {xj,tj) —)■ {x,t) as j ^ oo, tj > t and 

lim d{xj,tj) > d{x,t). 

j^ao 

We now argue in a similar way to the proof of (2), so that we obtain 05.21) for any z G Nt{x). 
Therefore 

OO 
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and it is easily seen that there exist 5 > 0 snch that 


B,{z) X + 5] C U 
t=i 

We thus conclude that z is an extinction point, and hence Nt{x) \Et = (/). □ 

Remark 5.1. (1) Theorem 15.41 (151) implies that if every x G T^ with t G (0,T) is a non 

extinction point then the distance function d is continuous in R” x (0, T) and hence 
Theorem 12.11 (Iml) holds. 

(2) If d is discontinuous at times 0 < ti < t 2 < ■■■ < < T (at one or more points in R"), 

we can apply Theorem 12.11 in the intervals (ti,t 2 ), ■■■, (tm,T). More precisely, 

under the assumptions of Theorem 12.11 we can show 

—)■ d locally uniformly in R"' x ((0, ti) U (fi, ^ 2 ) U ... U (tm, T)). 

9^+00 

Example 15.II in the next subsection shows that we can construct an evolution 

for which the associated distance function has discontinuities for each f G Q. Therefore, 

the idea described above cannot be applied. 

The next proposition gives a sufficient condition for the non-extinction condition. 

Proposition 5.5. Lett G (0,T). If x G int(rt), then x ^ Et. 

Proof. Let x G int(r 4 ). Then there exists a sequence C int(rt) that converges to x as 

j —)■ cx). Set Ej := dist(xj, dYf), which converges to 0 as j —)■ 0. Since we have B^.{xj) C T^, 
Proposition 15.21 (2) implies that C T. In particular Xj G Tt+(£j/L 2 )^ which is the vertex 

of the cone, and consequently we see that x is a non-extinction point. □ 

Remark 5.2. The converse of the assertion of Proposition 15.51 is not true in general. In 
fact, it is easy to construct the interface such that P* = {0} for all t G (0,T). Any x G Pj is 
a non-extinction point, but int(Pt) = 0. 

Remark 5.3. The opposite notion of an extinction point is an emerging point, which is 
dehned as follows: Let x G P*. We say that x is an emerging point if there exist e, 5 > 0 such 
that Bi;{x) X [t — S,t) C D. However, the property of hnite propagation implies that there 
are no emerging points. Suppose that x G P* is an emerging point, i.e., B^{x) x[t — S,t) <Z D 
for some e, 5 > 0. Choose M >1 large so that x G ^ subset of D 

by Proposition 15.21 ([ID. Thus x E D, a, contradiction. 

5.3 An Example 

In this subsection we present an example which shows that the idea presented in Remark 
15.11 d^D can not be applied, even if we restrict the evolutions to move inside a bounded domain 
instead of R". 
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Example 5.1 (A zero level set vanishing for all f G Qfl (0, T)). Let Q fl (0, T) = {ti, ^ 2 , •••}• 
Case 1. In R". 

Consider disjoint cubes with sides of length at least for n = 1, 2,.... Then inside every 
each one of them, we £t a circle Bt^ of radius The evolution of these circles under the 
equation 


(5.3) 


V = -l 


where V is the normal velocity (with normal pointing to the exterior of the circles), is given 
by 



Here R{t) is the radius of the circles. Notice that the evolution of (15.31) is the same as the 
zero level set of the solution u of the problem 


Ut = —|Vm| in R” X (0,T), 

u{x, 0) = uo(x) in R'^ 


if uq is for example the signed distance function to the circles Bn, with positive values in the 
interior of the circles. For a proof of equivalence of the two evolutions see for example m 
Section 4.2.3 and 4.2.4]. Then R(t) = tn — t and the extinction time of the circles is t = 

Case 2. In a bounded domain. 

Let H be a bounded open set. For every n G N we can hnd points G fl and positive 
numbers £„ such that B^^{xn) C f2 with B^^Xn) D B^^{xm) = 0 for n 7 ^ m. Then for 
On = £n/‘^iT + 1/2) we have := < enl‘2, and Bji^{xn) C B^{xn)- We then dehne 


Cn{x) 


Ctn 

O I rp rp 

c-n 

0 


in B,^/ 2 {Xn), 

in B,^{Xn)\B,^/ 2 {Xn), 
else 


and the velocity 

c(x) = sup Cn{x), 

n 

As in Case 1 we consider the problem 

Ut = —c(a:)|VM| 
u{x, 0) = Uq{x) 


for X G R”. 


in R" X (0, T), 

in R”. 


Here Uq is the signed distance function from the set UnGN with positive values in 

each Bji^{xn)- The extinction time of dBji^{xn) is as in the hrst case t = tn- 


6 Homogenization 


We conclude this paper by proving Theorem 12.31 Let us consider 


( 6 . 1 ) 


ut = Hi {X, 


1 + 9 


, Vm 
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and 


( 6 . 2 ) 


Ut = H2{u, Vm), 


where 6 = ^ 2/^1 is as in 02.201) . By the assnmptions on Hi and H 2 , the classical comparison 
and existence results still hold for the problems 06.Ih . 02 . 1 bl) and 06.21) . 02 . 1 bl) . 

To solve the problem 02 . 211 ) . 02.1bl) we use the notion of the iterative solution which was 
introduced in Remark 12.41 By the comparison and existence results for 06.11) and 06.21) . we 
see that 02 . 211 ) . 02 . 1 bj) admits a unique continuous iterative solution. 


6.1 Hamiltonians discontinuous in time 


Since the Hamiltonian H 12 is now discontinuous with respect to time, we have to be 
careful about the proof of our homogenization result. We do not use the notion of viscosity 
solutions introduced in Dehnition 12.41 where the upper- and lower semicontinuous envelopes 
are used for the equation, because otherwise we could not estimate the difference between 
{H 12 )* and (ili 2 )*. Thus we hrst discuss removability of the upper- and lower star of the 
equation as well as a connection between the iterative solution and the different notions of 
viscosity solutions of fl 2 . 2 ip . 

In this section we call u a viscosity subsolution (resp. supersolution) with star if it is 
a viscosity subsolution (resp. supersolution) in the sense of Dehnition 12.41 Also, we say 
that M is a viscosity subsolution (resp. supersolution) without star if it satishes the viscosity 
inequality 02.241) with F instead of F* (resp. F*). Note that, since F^ < F < F*, a 
viscosity subsolution (resp. supersolution) without star is always a viscosity subsolution 
(resp. supersolution) with star. Namely, a notion of viscosity solutions without star is 
stronger than that with star. 


Theorem 6.1. Let H be the iterative solution of 02 . 211 ) . 02.1bp . 
(1) H is a viscosity solution of 02 . 2 ip . 02 . 1 bP without star. 


(2) If V is a viscosity solution of 02.2ip . 02.1bp with star, then 


_ 


T)n 




Theorem 16.11 (fT]) asserts that H is a viscosity solution in R" x (0,T) not only in the 
sense with star but also in the sense without star. In other words, existence of solutions 
is established in both the cases. On the other hand, ([2]) is concerned with uniqueness of 
solutions since it asserts that any solution should be equal to u^. In the sense with star. 
Perron’s method (see Theorem 12. 6 p gives a viscosity solution up of fl 2 . 2 ip . fl2.1bp which is 
not necessarily continuous. By ([2]) we see that up = and therefore up is also a viscosity 
solution of fl 2 . 2 ip . fl 2 . 1 bp without star and an iterative solution as well. 


Proof. ([I]) We apply the fact that the viscosity property is extended up to the terminal 
time ([HI Section 7]). Since H is a viscosity solution of fl 6 .ip in R"" x (0, kiAt), we see that 
w^|R’*x(o,fciAt] is a viscosity subsolution of fl 6 .ll) in R"' x (0,fciAt]. This implies that H is a 
viscosity subsolution of fl 2 . 2 ip without star on R" x {kiAt}. Arguing in the same way on 
R" X {t} with t = e,e + kiAt, 2e,..., we conclude that H is a viscosity subsolution of (l2.21l) 
without star. The proof for supersolution is similar. 
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([2]) Since v and are, respectively, a viscosity subsolution and a supersolution of fl 6 .ll) in 
R" X (0, kiAt), the comparison principle for fl 6 .ip implies that v* < in R” x (0, kiAt). If 
we prove v* < on R” x {kiAt}, we then have v* < in R” x (fciAt, e) by the comparison 
principle for 06.21) . Iterating this argument, we hnally obtain v* < in R"' x (0,T). In the 
same manner, we derive < n* in R” x (0, T), and hence = v in R” X (0,T). 

It remains to prove that v*{x, kiAt) < u^{x,kiAt) for x G R"'. We now use the fact 

that u* is left accessible m Section 2, 9]), i.e., there exists a sequence such 

that tj < kiAt for all j > 1, {xj,tj) {x,kiAt) and v*{xj,tj) —)■ v*{x, kiAt) as j —)■ cx). 

Therefore, taking the limit in v*{xj,tj) < u^(xj,tj) gives v*{x, kiAt) < u^{x, kiAt). □ 

Remark 6.1. The same argument in the proof of ([2]) yields the comparison principle for 
fl2.2ip . Namely, if u and v are, respectively, a subsolution and a supersolution of fl2.2ip with 
star such that m*(-, 0) < u*(-,0) in R”, then u* < u* in R"" x [0,T). Also, similar arguments 
allow us to prove a local version of the comparison principle. Let (x, t) G R"' x (0, T) and 
r > 0. If u and v are a subsolution and a supersolution of fl2.2ip with star in Br{x) x {t — 
r,t + r) =: C, respectively, with u* < u* on dpC, then u* < u* in C. Here by dp we denote 
the parabolic boundary, that is, for H C R"' and a < b, 

dp{kl X (a, b)) := {dkl x [a, 6)) U (H x {a}). 

Remark 6.2. See |22l [5] for more results concerning Hamiltonians discontinuous in time. 

6.2 Cell problems 

We study an one-dimensional cell problem with discontinuity, whose solution and eigen¬ 
value will be needed in the proof of our homogenization result. Consider 

(6.3) v'{t) + X = H{t) in T, 

where T = R/Z is the one-dimensional torus, H G T^(T) and A G R. Although we only 
need to study piecewise continuous H for our homogenization result, we here take it as a 
L^-function since the technical aspects of the proof allow us to generalize H without any 
additional effort. For the special case where H is piecewise continuous, see Remark 16.31 We 
dehne 

:= limsup [ H(s)ds) , H#(r) := limini (y [ H(s)ds) . 

^4^0 Vfc Jr—k J \k Jr—k / 

Lemma 6.2 (Solvability of the cell problem). We set 

(6.4) A := / H{s)ds, v{t) := v{0) — Xr + [ H{s)ds. 

Jo Jo 

Then v is a viscosity solution of 06.31) in the following sense: //maxT(u — (f) = {v — 0)(ro) 
(resp. minT(u — (f) = {v — 0)(ro )) for tq G T and cf) G C^(T), then 

(6.5) (j)'{To) + X< H#{to) (resp. 0'(ro) A > H*{to)). 
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Proof. 1. We first note that n is a periodic function thanks to the choice of A. Indeed, for 
all r G R and m G Z, we observe 


v{t + m) 


pT-\-7n 

n(0) —A(r + m)+ / H{s)ds 
Jo 

v{0) — + m) + (^Xm + J H{s)ds^ 

n(0) —Ar+ / H{s)ds 
Jo 

v{t). 


Thus V is periodic. 

2. Take tq G T and 0 G C^(T) such that maxT(n — 0) = (w — 0)(ro). For fc > 0 we have 


0 (^o) - <Piro -k) . v{to) - v{to - k) _ 1 

k - k -~^^k 


rro 


'ro—/c 


H{s)ds. 


Taking hminffc|o implies the hrst inequality in (16.51) . A similar argument shows that n is a 
super solution. □ 

Remark 6.3. Let 0 = Tq < Ti < • • • < tat = 1 be a partition of [0,1] and assume that 
H G T^(T) is continuous on each (tj, Tj+i]. Then we have = H, and consequently 

the viscosity inequalities in 06.51) become 

(p'{ro) + A < H{to) (resp. 0'(ro) + A > H{to)). 

In other words, v given by 06.41) is a viscosity solution of 06.31) without star. 


6.3 Proof of homogenization 

Proof of Theorem \2.,J\. 1. Let be the iterative solution of 02.211) and 02.1bl) . We denote 
by u and u the upper- and lower half-relaxed limit of respectively, i.e., u = limsup*^QM^ 
and u = lim inf*£_,.o Since the functions Uq{x) — Kt and Mo(a;) + Kt with K > 0 large 
are, respectively, a subsolution and a supersolution of 02 . 2 ip . it follows from comparison 
that Uo{x) — Kt < K{x,t) < Uq{x) + Kt. This implies —cx) < u < u < -|-cx) and 0) = 
u{x, 0 ) = uo{x). 

2. Let us show that m is a subsolution of 02.22p . Let 0 be a test function for u at {x,i) 
from above, i.e, 

(6.6) u < (f) in {Bji{x) X {i — R,i + R)) \ {{x,i)}, 

(6.7) u{x,i) = (p{x,f) 

for some R > 0 such that 0<f—i?<f-|-i?<T. We set 

H{t) ;=Ri2(£,t>,0,V0), 
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where 0 = 0(f, t) and V0 = V0(f, t). Then H G -h^(T) and H is continnous on (0, ki^t/e] 
and (fciAt/e, 1], Let v and A be as in fl6.4p . By Remark 16.31 we see that n is a viscosity 
solntion of (Ib.dl) withont star. Noting that 6 = ^ 2 /^ 1 , we observe 


Hu{x,t,T,r,p)dT = I " Hi\x,^—^,p] dr + j H 2 ir,p) dr 


This implies that 


ki 


ki + ^2 
1 


Hi \x,Y^^p\ + 


1 + 9 

= H{x,t,r,p) 


ki + k2 

(^’ r^’^) ^ ^H2{r,p)j 


H2ir,p) 


A = / H{T)dT= / Tri 2 (x, t, r, 0, V0) dr = iL(x, t, 0, V0). 


Consequently, v solves 


(6.8) n'(r)+id(f, t, 0, V0) = idi 2 (r,f, r, 0, V0) in T. 

3. We want to show that 

0f < H{x, i, 0, V0) 

with 0t = (j)t{x, t). Suppose in the contrary that there is /i > 0 such that 


(6.9) 


0t > id(f,t, 0, V0) +/i. 


Let us introduce a perturbed test function. Dehne 

(j)^{x, t) = 0(r, t) + ev . 


Since v is bounded, we see that 0^ converges to 0 uniformly as £ —?■ 0. We will show that 0^ 
is a supersolution of fl2.2ip in Br{x) x (t — r,t + r) =: C, where r G (0, R) is chosen to be 
small so that 


( 6 . 10 ) 

( 6 . 11 ) 


Hi 2 ( X, t, (j){x, t), V0(a:, t) 


-H 


12 


x,i, j,0, V0 


t\ - 


< 




for all {x,t) G C. Although H 12 = Hi 2 {x,t,T,r,p) is discontinuous in r, fl6.1ip is achieved 
because we £x r = to/e. More precisely, 06.111) is satished if 


Hi (x, t, V0(x, t)) — Hi{x, i, V0) 



H 2 (0(x, t), V0(x, t)) - iL2(0, V0) 


< 


p 

8 
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for all (x, t) e C. Allowing a larger error, we are able to replace 0(x, t) on the left-hand side 
of fib.lip by (j)^{x,t) with £ > 0 small enough. Namely, we have 


( 6 . 12 ) 


Hi2 (x,t, —,0^(x,t), V0(x,t) ) - Hi2 (x,t, —,0, V0 




4. Let -0 be a test function for 0^ at (xo,to) ^ C from below. Then the function 


r v{t) - - er) - 0 (xo, er)) 

has a local minimum at Tq := to/e. Also, from the smoothness of it follows that 


(6.13) 


V0^(xo,to) = V(/)(xo,to) = V?/>(xo,to)- 


Since f is a viscosity super solution of (l 6 . 8 l) . we have 

— ^ /V /V f 

'ipt{xo,to) - (j)t{xo,to) -h i^(x,t,0, V0) > Hi 2 [x,t, —,(j),V(j)] . 

Let (j)l = (j)^{xo,to), V0O = V0(xo,to) and V^/’o = V?^(xo,to)- Applying (16.101) . (16.121) and 
(16.op to the above inequality, we compute 


(6.14) 


Aixo, to) > 0t(xo, to) - H{x, t, 0, V0) -h Hi 2 i, j, 0, V0^ 

> 0t - ^ - H{x, t, 0, V0) + Hi 2 (^Xo, to, 00, V0O^ - ^ 

> Hi 2 ^Xq, to, —, 00 V 0 o^ + 


For the last inequality we have used 06.131) . The above inequality shows that 0^ is a super¬ 
solution. Moreover, since H 12 = Hi 2 {x,t,T,r,p) is continuous in the r-variable, the estimate 
06.14P implies that there is a small r/o > 0 such that 0^ — 77 is also a supersolution of 02.211) 
for every 7] G ( 0 , rjo]. 

5. Set 


5o := -max(M - 0), 

opC 


which is positive by 06.6p . Also, let 6 := min{(5o/2, rjo}. We then have 


max(M^ — 0^) < —6, 

dpC ^ -r j - 

i.e., < 0^ — (5 on dpC for e > 0 small enough. We now apply the comparison principle 

for a subsolution and a supersolution 0^ — 5 of 02.2ip to obtain < 0^ — 5 in C. Taking 
limsupg^Q at (x,t), we see u{x,i) < 0(x, t) — 5. This is a contradiction to 06.7p . and hence 
M is a subsolution of 02 .22p . 

6 . Similarly we show that m is a supersolution of 02.221) . and therefore u = uhj compar¬ 
ison. This implies the locally uniform convergence of to the unique viscosity solution 
of arm and fl2TB . □ 

Remark 6.4. As long as the comparison principle is true, this homogenization result still 
holds for more general equations with Hi and H 2 which are not necessarily of the forms 
Hi = Hi{x,t,p) and H 2 = H 2 {r,p). 


39 































A Lipschitz continuity of solutions 


The properties of the solution of the problem (12.61) and fl2.1bll might come in handy 
when studying numerical results. For this reason we prove here a Lipschitz estimate for the 
solution, under the assumption that the initial datum is Lipschitz continuous. We also give 
an explicit representation of the Lipschitz constant in terms of the Lipschitz constant of the 
initial datum and the Lipschitz constant of the Hamiltonian Hi denoted by D[t) as in flH3-sj) . 
Although there are plenty of results in the literature concerning the Lipschitz continuity of 
viscosity solutions, a Lipschitz estimate for the Hamiltonians which are being studied in this 
paper does not exist up to the authors’ knowledge. Moreover we are more concerned in a 
Lipschitz constant that does not depend on the parameter 6 . 

Proof of Pro‘position \2. 71 1. Let $(a:, y, t) = u{x, t) — u{y, t) — L(t)\x — y\ for x,y E R" and 
t G [0,T). We proceed by contradiction. Suppose that 

M= sup ^{x,y,t)>0. 

x,yeR",iG[0,T) 

Since u has at most linear growth (Theorem 12.61) and uo is Lipschitz continuous with 
Lipschitz constant Lq, we have 

u{x, t) - u{y, t) < Uo{x) +Kt- {uo{y) - Kt) 

(A.l) <CT + LQ\x-y\ 

with Ct = 2KT. We dehne 

^a{x,y,t) = u{x,t) - u{y,t) - L^{t)\x - y\ - - a{\x\‘^ + \y\'^), 


where 


Set 


L„(f) = max |Lo, 1 + 2 ^^'^ g/o 
p(t) = 2\JaCTD(t). 

= sup ^a{x,y,t) 

a:,yGR"pG[0,T) 


for a = {t], a). Since u has at most linear growth, there are x^, ya G R" and G [0, T) such 
that 

(^CT ) 2/t7) ^(t) • 

2. By the dehnition of M, for every h > 0 there are xs, ys and ts such that 


^{xs,y5,t5) > M - 5. 

C t 

Since La{t) —?■ L{t) = max{Lo, Ije^o uniformly in f as a —?■ 0, there is £ > 0 small 

enough, independent of h, such that 

-L^{ts) > -e- L{ts). 
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Then for 5 = M/4 we have 

Ma > ^a{xs,y5,ts) 

Tj 

= <h(a; 5 , ys, ts) - -—- - a{\xs\^ + lys]^) - £\xs - ys] 

I — Is 

> M -6 - ^ - a{\xs\^ + \ys\^) - e\xs - ys\ 

I — ts 

M 

>y>0 

for T], a, e small enough, since 6 is fixed. From this it follows that 

(A.2) ^^{xa,ya,U) > 0. 

3. We claim 

(A.3) a\xa\, a\y^\ < \[aC^. 

By flA.ll) we observe 

u{x,t) - u{y,t) - La,{t)\x -y\<CT+ (To - Lo,{t))\x - y\ < Ct 
for all (x, y, t), and hence we can write 

u{x„,U) - u{y„,U) - La{t„)\x„ - y„\ - —-< Ct- 

1 tfj 

The left-hand side is equal to ^^{xo-,ya,ta) + a\xa\^. By (IA.2I1 we get 

a|xo-|^ < Ct- 

Similarly we have < Ct, and these inequalities show flA.3F 

4. We prove > 0 and x^ 7 ^ ya- Suppose that to- = 0. Then, since uq = m(-, 0) is 
Lipschitz continuous with Lipschitz constant Lq < Tq,( 0), we have 

0 < 1 /^, 0 ) < u{x„,0) - u{y^,0) - T„( 0 )|a:<, - y^\ < 0 , 

a contradiction. Since ^a{xa,Xa,ta) < 0, we have that x^ 7 ^ y^j- 

5. Since x^ 7 ^ y^, we have |x — |/| > 0 in a neighbourhood of {xa,ya)- Therefore, we can 
apply HB Lemma 2] and get, for = {x^r — ya-)/\xa- — y^l, 

^ + L'a{ta)\Xa - ya\ 

C {Hi(^Xq^, f(7, -\- ‘2(y.X(j^ f^-, LQi(tu')Pa ^crpo-)} 

-h 9{(3{u{x^,Oh{La{C)p^ -h 2 q!X<^) - (3{u{y^,0)h{LaiC)p^ - 2ay^)} 

(A.4) =: Ji -|- 12- 

We can also rewrite I 2 as 

h = 0 {(3{u{x^,0) - (3{u{y^,0)} h{La,{C)p^ + 2ax^) 

(A.5) 9/3{u{ya,C)) ■ {h{La(t^)p^ + 2ax„) - h{La{C)p„ - 2ay„)]. 
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6. Let us give estimates of Ji and 12- Using flH3-sD . the Lipschitz continuity in x of iLi, 
together with flH2p . we get 

h < D{t^)\x^ - y„\ ■ \La{t„)p^ + 2ax^j 

-j- , tfj-y 2crXo-) Lo^ito-^Pf^ 2o|/o-)}- 

The hrst term on the right-hand side can be estimated by (IA.3P and |po-| = 1, while we apply 
(IH4H . the Lipschitz continuity in p of iLi, to the second term. Then 

h < D{U)\xa - Hal {La{t„) + 2^aCT^ + 2aL2\xa + ya\- 
By the dehnition of p we have 

(A.6) Ji < \x^ - y„\ {D{t^)La{ta) + piU)) + 2aL2\x^ + y^\. 

We next show that the hrst term on the right-hand side of flA.5l) is not positive. We hrst 
note that u{xaito) > u{ya,ta) by flA.2l) . Therefore, we have (3{u{xa,ta)) > (3{u{ya,ta)) since 
fi is increasing. As a second remark, by the dehnition of La and flA.3ll . we have 

\La{ta)pa + 2Q!a;^| > - 2 ^/aC^ > 1. 

Then the assumption (12.101) and the continuity of h imply that h{La(ta)pa + 2aXa) < 0. 
According to these two remarks we have 

h < 9\(3{u{ya,ta))\ ■ \h{La{ta)Pa + 2aXa) - h{La{ta)Pa “ 2a|/^)|. 

Using the boundedness of (3 and the uniform continuity of h, 02.81) . we can further estimate 
the above as follows: 

(A.7) l2<0Muh{2a\xa + ya\), 

where M is an upper bound for |/9|. 

7. Applying 0A.6D and 0A.7D to 0A.4D . we get 

(A.8) ^ + L'^{ta)\Xa-ya\ < \Xa - ya\{La{ta) D{ta) + pit^)) + J, 

where 

J = 2aL2\xa + ya\ + 9Mu:h{2a\xa + p^l). 

Note that the function La has been chosen so that it solves the diherential equation 

According to this, the estimate 0A.8P becomes 

^ < h(^<7)(l - La{ta))\Xa - | J. 

Since Lq, > 1, we have 

^ < J = 2aL2\xa + ya\ + 9Mujh{2a\xa + y^D- 
Using 0A.3D . we can send a ^ 0 and get a contradiction. □ 
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